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BOUNDS  ON  DOPPLER  FREQUENCY  ESTIMATION  IN  CORRELATED 
NONSTATIONARY  GAUSSIAN  NOISE 


INTRODUCTION 

There  exist  many  radar  problems  that  require  the  detection  of  a  signal  in  additive  Gaussian 
noise.  The  optimal  detector  for  this  case  is  the  classical  matched  filter,  which  is  obtained  by  applying 
the  Neyman- Pearson  criteria  to  a  binary  hypothesis  test  [1].  For  the  pulsed  radar  problem  of  detect¬ 
ing  a  single  Doppler-shifted  target  echo  in  a  given  range  cell,  there  are  three  unknowns  in  the 
matched  filter:  the  covariance  matrix,  the  complex  signal  amplitude,  and  the  target  Doppler  shift 
present  in  the  steering  vector.  One  practical  detector  is  usually  implemented  as  follows:  the  covari¬ 
ance  matrix  is  estimated  from  data  in  nearby  range  cells  [2];  a  bank  of  Doppler  filters  or  steering  vec¬ 
tors  is  implemented  so  that  the  signal  falls  in  one  of  them;  finally,  the  magnitudes  of  all  the  Doppler 
filters  are  found  and  compared  to  a  threshold.  This  final  step  eliminates  the  need  to  know  the  signal 
phase,  and  the  signal  magnitude  is  not  required  because  the  resulting  test  is  uniformly  most  powerful. 
Reference  3  describes  a  variation  of  this  procedure. 

After  detection,  the  Doppler  shift  of  the  target  and  its  complex  magnitude  is  estimated  Refer¬ 
ence  4  gives  the  maximum  likelihood  estimate  (MLE)  of  the  complex  magnitude,  and  the  present 
report  develops  the  estimation  of  the  Doppler.  The  performance  of  the  MLE  of  Doppler  is  studied 
for  representative  cases  of  interest  that  contain  Gaussian  noise  composed  of  thermal  noise  plus  clutter 
plus  uncorrelated  interference  and  a  Doppler-shifted  single  point  target.  The  performance  of  the 
MLE  is  then  compared  with  that  of  the  Cramer-Rao  bound. 

MAXIMUM  LIKELIHOOD  ESTIMATE 

The  received  signal  vector  X,  composed  of  n  complex  samples,  is  given  by 

X  =  N  +  bW , 

where  N  is  a  complex  zero-mean  Gaussian  process  with  covariance  R,  b  is  the  unknown  complex 
signal  strength,  and  W  is  the  steering  vector.  The  steering  vector  is  written  in  terms  of  A,  the  phase 
shift  from  pulse  to  pulse,  owing  to  Doppler  as 

WT  =  [1  ejA  . . .  eJ{n  ~  nA], 

where  T  denotes  the  transpose,  and  j  is  the  square  root  of  -  1 .  The  joint  probability  density  of  X  is 
then 

p(X|0)  =  exp  \-(X  -  bW)T  R~\X  -  bW)\ ,  (1) 

ir  \  R\ 

where  0  represents  the  unknown  parameters  A,  b,  and  b ,  and  (•)  denotes  the  complex  conjugate. 
The  MLEs  of  A,  b ,  and  b  are  obtained  by  differentiating  the  likelihood  density  with  respect  to  each 
of  these  variables,  setting  these  resulting  equations  equal  to  zero,  and  solving  for  the  three  variables. 
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The  MLEs  of  b  and  b  are  well  known  [4]  and  are  given  by 


WtR-{X 
WtR-'W 
and 

r  =  XTR~lW 

wtr-'w' 

The  partial  derivative  of  Eq.  (1)  with  respect  to  A  yields 

r  aw  I 

+ 


.  dW 
* 


(X  -  bW)T  R~x 
Substituting  Eqs.  (2)  and  (3)  into  Eq.  (4)  yields 

(WTR~lW)  (WTR~lX) 


r  -r 

r  dWT 
b 


dA 


R~\X  -  bW)  =  0. 


.-To-,  aw 


W'R 


dA 


(2) 


(3) 


(4) 


+  (WTR~lW) 


aw‘ 

dA 


RlX 


(XrR~lW) 


WTR- 

,  dW_ 
dA 

dWT 

...1 

(WtR-'X)  ( X‘R-'W ) 


TD-h 


dA 


(wtr~'X)  ( xtr~'w )  =  0. 


which  can  be  reduced  by  using  the  definition  of  W  to 

£  (q  ~  k  +  i  -  /)  R« XmXp  =  0.  (5) 

i  ,k  ,1  ,m  ,p  .q  =  1 

The  root  of  this  polynomial  that  maximizes  p(X  |  0)  in  Eq.  (1)  is  the  MLE  of  A. 

CRAMER-RAO  BOUND 

The  MLE  is  now  compared  to  the  Cramer-Rao  bound  [5,6].  The  bound  on  the  variance  of  A  is 
obtained  from  the  Fisher’s  information  matrix  J  whose  elements  are 


4  =  £ 


/• 

r  ^ 

ainp(X|0) 

d\np(X\G) 

> 

ae, 

l  J 

(6) 
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where  dx  =  A,  02  =  b,  03  =  b;  E  is  the  expected  value  operator;  and  i  and  j  are  the  indices. 
By  using  the  density  function  (Eq.  (1))  and  the  definition  (Eq.  (6)),  then 

,  awT  dw  ,  L-  dwT  dw 

11  dA  dA  dA  dA 

J2  2  =  WtR~1W  =  C, 


y33  =  wtr~1w  =  c. 


*  * 


a  vltT  _ 

Jn  =  b  ~  R-1  W  =  bA , 


—  dW*  _ 

7 i3  =  b  ^-R-'W  =  bA, 


dA 


and 


72i  =  bWTR~l  ^  =  b  =  bA, 

dA  a  A 

731  =  bWTR~l  ~  =  b  /?"'  iy  =  bA, 
31  dA  dA 


•^23  ~  -7  32  ~  0. 

The  matrix  7  can  then  be  inverted  in  a  closed  form 


C  C 

-bAC 

-bAC 

7-1  _  1 

-bAC 

CJU  -  bbAA 

b2A2 

|7  | 

-bAC 

Pa2 

CJ  ] )  -  bbAA 

(7) 


where  1 7  | ,  the  determinant  of  the  matrix  7 ,  is  expressed  by 

|7  |  =  CCJU  -  bbAAC  -  bbAAC . 

The  bound  on  the  variance  of  the  estimate  of  A  required  by  the  Cramer-Rao  bound  is 

2  CC 
—  TTT  ’ 


(8) 


which  is  the  upper  left-hand  corner  of  the  matrix  given  in  Eq.  (7). 
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A  RADAR  ESTIMATION  PROBLEM 

The  covariance  matrix  R,  composed  of  the  sum  of  three  covariance  matrices,  is  written  as 

R  —  Rg  Rfl  -E  Rj  , 

where  Rc,  RN,  and  Rf  are  the  covariance  matrices  of  clutter,  thermal  noise,  and  interference,  respec¬ 
tively.  For  the  examples  to  be  studied  these  are  given  by 


1 

p 
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On 

where  ac,  aN,  and  a,s  are  the  standard  deviations  of  the  clutter,  thermal  noise,  and  the  ith  interfer¬ 
ence  component,  respectively,  and  p  is  the  Markov  correlation  coefficient.  The  signal-to-noise  ratio 
(SNR)  is  defined  as 


SNR  =  4- 

oN 


The  performance  of  the  MLE  is  evaluated  by  using  simulation.  Random  numbers  are  generated 
according  to  the  equation 


X  =  Nn  +  Nc  +  N,  +  |  6  |  exp  (;<*>)  A=Ao, 

where  fV  is  evaluated  at  A  =  Ay,  the  true  Doppler  shift;  Nn,  Nt ,  and  Nc  are  Gaussian-distributed 
complex  vectors;  and  0  is  a  random  phase  distributed  uniformly  between  0  and  2ir.  The  correlated 
clutter  samples  are  generated  from  the  uncorrelated  Gaussian  random  numbers  given  by  Z  using 
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where  Rc  =  LLT  and  L  is  a  lower  triangular  matrix.  For  the  3  x  3  cases, 


These  random  numbers  are  used  in  Eq.  (5)  to  obtain  the  MLE  of  A  per  trial.  The  variance  of  these 

naA 

estimates  is  then  calculated  by  using  a  number  of  trials.  The  variances  are  normalized  by  — —  and 

2ic 

are  plotted  vs  SNR  in  various  clutter  and  interference  environments. 

The  Cramer-Rao  bounds  of  the  normalized  variances  can  also  be  obtained  as  a  function  of  SNR 
for  various  cases.  Equation  (8)  shows  that  the  logarithm  of  the  normalized  Cramer-Rao  bound  is  a 
linear  function  of  SNR  in  dB. 

EXAMPLE  RESULTS 

The  results  of  the  MLEs  for  the  case  of  n  =  3,  p  =  0.999,  and  a^/a^  -  20  dB  without 
interference  are  shown  in  Figs.  1  and  2  (in  dashed  lines)  for  A  =  O.lir  and  A  =  0.5ir,  respectively. 
These  results  are  very  close  to  the  Cramer-Rao  bounds  (the  solid  lines  in  the  figures)  when  the  SNR 
is  comparable  to  or  above  the  clutter  strength.  Figures  3  through  7  show  the  results  for  the  case 
where  n  =  3,  p  =  0.999,  oc/°n  =  20  dB,  and  A  =  0.5ir  in  the  presence  of  interference.  Again, 
the  MLE  errors  are  tightly  bounded  by  the  Cramer-Rao  bounds  for  SNR  equal  to  or  greater  than  the 
stronger  of  the  clutter  and  the  interference  when  only  one  pulse  is  subject  to  interference  (Figs.  3  and 
4).  As  the  interference  is  applied  to  more  than  one  pulse,  the  required  SNR  for  the  MLE  errors  to  be 
close  to  the  Cramer-Rao  bounds  is,  generally,  higher  (Figs.  5  to  7). 

Figure  8  shows  the  results  for  the  case  of  n  =  5,  p  =  0.999,  and  A  =  0.5ir  without  interfer¬ 
ence  in  the  presence  of  clutter.  The  MLE  errors  are  nearly  identical  to  the  Cramer-Rao  bounds.  Fig¬ 
ure  9  shows  the  results  for  a  clutter  situation  that  has  a  clutter-to-noise  ratio  of  20  dB  when  interfer¬ 
ence  is  applied  to  both  the  first  and  the  second  pulses.  Again  the  two  methods  yield  very  similar 
results. 

SUMMARY 

After  detecting  a  radar  target  by  use  of  a  matched  filter  in  the  presence  of  clutter,  thermal  noise, 
and  uncorrelated  interference,  the  target  Doppler  shift  is  estimated.  The  MLE  of  Doppler  and  the 
associated  Cramer-Rao  bound  were  found  for  the  Gaussian  noise  case.  The  estimation  performance 
(variance  of  the  estimate)  was  found  for  several  cases  of  interest.  The  estimates  were  typically  near 
the  Cramer-Rao  bound.  In  fact,  they  could  be  quite  close  even  in  the  presence  of  high  clutter  and 
interference. 
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CRAMMER-RAO  BOUND  - 

MAXIMUM  LIKELIHOOD  ESTIMATE - 


Fig.  1  —  The  bounds  and  the  maximum  likelihood  estimates  of  normalized 
Doppler  shift  errors  vs  SNR  for  n  =  3,  p  =  0.999.  and  A  =  0. 1  tt 
without  interference  in  clutter 
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CRAMMER-RAO  BOUND  - 

MAXIMUM  LIKELIHOOD  ESTIMATE - 


Fig.  2  —  The  bounds  and  the  maximum  likelihood  estimates  of  normalized 
Doppler  shift  errors  vs  SNR  for  n  =  3,  p  =  0.999.  and  A  =  0.5  rr 
without  interference  in  clutter 
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CRAMMER-RAO  BOUND  - 

MAXIMUM  LIKELIHOOD  ESTIMATE - 


Fig  5  —  The  bounds  and  the  maximum  likelihood  estimates  of  normalized 
Doppler  shift  errors  vs  SNR  for  n  =3,  p  =  0.999.  A  =  0.5  ir.  and 
ol/ol  =  20  dB  with  interference  on  the  first  and  the  second  pulses 
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CRAMMER-RAO  BOUND  - 

MAXIMUM  LIKELIHOOD  ESTIMATE - 


SNR  (dB) 

Fig.  6  —  The  bounds  and  the  maximum  likelihood  estimates  of  normalized 
Doppler  shift  errors  vs  SNR  for  n  =  3,  p  =  0.999,  A  =  0.5  t.  and 
Of /ct.v  =  20  dB  with  interference  on  the  first  and  the  third  pulses 
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SNR  (dB) 
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